Electronic dielectric constants of insulators by the polarization method 
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We discuss a non-perturbative, technically straightforward, easy-to-use, and computationally 
affordable method, based on polarization theory, for the calculation of the electronic dielectric 
constant of insulating solids at the first principles level. We apply the method to GaAs, AlAs, InN, 
SiC, ZnO, GaN, AIN, BeO, LiF, PbTiOa, and CaTiOa. The predicted e°°'s agree weU with those 
given by Density Functional Perturbation Theory (the reference theoretical treatment), and they 
are generally within less than 10 % of experiment. 
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The electronic dielectric constant £°°, measuring the 
response to a uniform electrostatic field, is a fundamental 
quantity in basic and applied solid state physics. Besides 
its intrinsic interest, knowledge of £°° is crucial to the 
calculation of the static dielectric constant e°. Ab initio 
calculations oie°° have been performed in recent jtimes by 
density functional perturbation theory (DFPT)trEl for a 
host of different materials. Here we discuss a new method 
to compute e°° based on density functional jtheoryo and 
the modern theory of dielectric polarization,El and apply 
it to a set of polar crystals. Despite its simplicity and ease 
of use, the method predicts e°°'s in close agreement with 
those obtained by DFPT. This suggests that theoretical 
predictions deviate from (and typically overestimate) the 
experimental value on account of density functional the- 
ory, and not of the specific method or implementation. 

We have sketched the basics of the method in a report 
on the calculation of the static dielectric constant B Al- 
most all the ingredients needed to evaluate the latter can 
be obtained from ab-initio calculations of total energy, 
stress, forces, and dielectriC|-polarization in zero field (by 
the Berry-phase approachaQ) for bulk systems. The no- 
table exception is the electronic dielectric constant 
for which we follow the alternate approach described be- 
low. 

Theory - The core of the argument is that £°° 
can be obtained from the relationship between macro- 
scopic polaripi^ipn in zero field and interface charge 
accumulatiorJja'tl in appropriate superlattices. An insu- 
lating superlattice is constructed, which consists of peri- 
odically alternating slabs of equal length, stacked along 
some fixed direction and made of materials 1 and 2. In 
such a superlattice, in the absence of external sources of 
fields, the displacement field orthogonal to the interfaces 
is conserved: 
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El + AttPi{Ei) ^E2+ 4nP2{E2) = D2. (1) 

We expand the polarization to first order in the screened 
fields in the two materials, indexed by i: 

P,(E) ^ Pf^ + X^E^, (2) 

where Pf is the zero field polarization and Xi fl^6 sus- 
ceptibility of material i. The presence of a zero field 



polarization is important: substituting the last relation 
into Eq. one obtains 
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(In this context, it does not matter whether the zero field 
polarization is spontaneous or piezoelectric, or a combi- 
nation of the two.) Only for a null P^ does one recover 
the familiar equality 



e^Ei = e^E2 . 



(4) 



To proceed further, one notes that periodic boundary 
conditions imply 

E ~ El — — E2-, /\E = El — E2 = 2ii/, 

and therefore Eq. || becomes 



4^(pf -p(")) = -(e-+e-)Ai?. 



(5) 



Recalling that the charge accumulation per unit area 
at the interface between materials 1 and 2 is Sint — 
±AE/4tt, we finally obtain 

Sint = ±2{P^"^ -P^°^)/{er+eT) 

= ±2 Ap("V(£r + (6) 

This relation connects the difference in macroscopic bulk 
polarization at zero field with the components ej°2 of the 
dielectric tensors-of the interfaced materials along the 
interface normal.ll3 

In an undistorted homojunction, i.e. a superlattice in 
which material 1 is identical to material 2, there is ef- 
fectively no interface. Therefore there is no polarization 
change, and the interface charge is zero. It is neverthe- 
less possible to generate a polarization difference in a 
controlled manner, by inducing a small distortion 5 of 
one of the atomic sublattices in half of the superlattice 
unit cell. This produces a difference in polarization, and 
a charge accumulation at the interface. The interface 
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charge Sint at the interface between distorted and undis- 
torted reapHS can be easily calculated via macroscopic 
averagingS'tll of the charge density. On the other hand, 
the zero field polarizations P2 for the bulk material in 
the undistorted state, and Pi for the bulk material in 
the same strain state as in the superlattice, are evalu- 
ated using the Berry phase technique. From Eq. (||), one 
then extracts the average electronic dielectric constant 
= {ef +e^)/2. 
By construction £2 equals the dielectric constant, while 
ei, the dielectric constant in the distorted state, of course 
does not. Therefore e does not equal the sought-after di- 
electric constant for (5 7^ 0; However, e does equals the 
dielectric constant (the tensor component along n) in the 
limit of zero distortion: 

e°° = Hm e. 

This limit can be evaluated with essentially arbitrary 
accuracy by extrapolation or interpolation. Also, since 
the materials involved in the heterojunction are identi- 
cal modulo a vanishingly small distortion, it is virtually 
guaranteed that no interface state exists, so that no band 
bending or electrostatic perturbation adulterates the po- 
larization effects. 

Summarizing, in the present approach the dielectric 
constant is simply obtained using the geometric quan- 
tum phase polarization and relatively small, accurately 
controllable supercell calculations. In the latter calcula- 
tions, the slabs should be short enough that the constant 
electric field will not cause metallization, and the slabs 
should be sufficiently long to recover bulk-like behavior 
away from the interfaces. Both requirements are gen- 
erally met also by materials with small calculated gaps 
if sufficiently small strains are applied. Unlike DFPT, 
the present method does not resort to perturbation the- 
ory, and it is novel in that the determination of the elec- 
tronic screening uses the connection with the geomet- 
ric quantum phase. E3 The calculations involved are non- 
intensive, and can be performed even on Pentium-like 
personal computers. Also, since the implementation is 
very much simpler than that of DFPT, the method seems 
promising as a general-purpose tool for non-specialists. 

Applications - We now apply the method to a set of 
representative polar materials of general interest: SiC, 
GaAs, AlAs, InN, GaN, AIN, BeO, ZnO, LiF, cubic 
PbTiOa, and cubic CaTiOa. All our calculations, as de- 
tailed below, are done in the local density approximation 
to density functional theory.a Our results are compared 
to results of DFPT calculations (where available) , which 
effectively function as reference for new computational 
methods in this area. Comparison with experiment is 
also presented when possible. 

In the calculations we use|— ■ the Ceperley- 
Alder exchang&a:orrelation energyll3 and ultrasoft 
pseudopotentialsO for the electron-ion interaction. The 
pseudopotentials have been constructed to include the 
following semicore states in the valence manifold: Zn and 



Ga 3(i, In 4(i; Li Is; Pb 5d, Ti and Ca 3s3p. A plane- 
wave basis with 25 Ryd cutoff is found to be sufficient 
to converge the quantities of interest in all the materials 
investigated. Bulk Brillouin zpae summations are done 
over appropriate Chadi-Cohenlj k-point meshes for the 
relevant structures. Bulk polarizations-are obtained in all 
cases via the Bepy phase techniqueEnZl using a 16-point 
Monkhorst-PackEEl k-point mesh in the a-plane direction 
and 10 point uniform mesh in the c direction. For the 
supercell calculation, we have employed superlattices in- 
cluding typically 16-20 atoms, oriented along (0001) for 
wurtzites, (111) or (100) for zincblende, (100) for the 
NaCl structure, and (100) for cubic perovskites. (Note 
that wurtzite has two independent components of the di- 
electric tensor, and the one we are actually calculating is 
that along the polar axis.) In the supercell calculations, 
downfolded meshes were used comprising 12 k-points for 
wurtzite and zincblende, 8 points for the NaCl-structure, 
and 12 points for perovskites. The cation sublattice dis- 
placements 5 are typically 0.05-0.1 % of the bond length. 
Ionic relaxation is never allowed, so that the response 
is purely electronic. All the calculations are performed 
at the theoretical lattice constants, that are reported in 
Table |. 

With the above reported theoretical ingredients, we 
obtained the theoretical values of e°° listed in Table || 
together with DFPT and experimental values, and plot- 
ted in Fig. versus the experimental values for the dif- 
ferent materials. The general level of agreement seems 
quite good on the scale of Fig. |^. To give a closer view, 
in Fig. H we display the relative percental deviation of 
the theoretical with respect to experiment, both for 
our method and DFPT. 

Discussion - The main content of Fig. || is that 
DFPT and the present method agree quite closely (the 
deviation for GaN is probably due to a different treat- 
ment of the Ga id electrons). It thus appears that de- 
viations from experiment are not related to the specific 
method used, but are likely to be a token of the underly- 
ing density functional formalism. Most theoretical values 
are overestimates of the experimental data, the main ex- 
ception, being the 4 % underestimate for PbTiOs, both in 
DFPTU (at the experimental lattice constant) and in the 
present method (at the theoretical lattice constant). This 
is possibly due to the uncertainties in the experim(pi|tal 
values, which are in fact plasmon-pole extrapolationscil to 
optical a;=0 of values measured in the visible at relatively 
high temperature, while the calculation is at zero tem- 
perature and zero lo (in fact, this holds also for CaTiOs). 
Indeed the situation for perovskites (even in the para- 
clectric phase) is far from settled in general, ather re- 
cent linear response results of another groupc£lEj over- 
estimate experiment considerably : the e°° reported fpt 
cubic SrTiOa is 6.63 compared with 4.69 experimental,E£l 
and also_6..63 compared with 5.18 experimental for cubic 
KNbOa.m 

It should be mentioned that £°° in our scheme is actu- 
ally a finite-g value due to the finiteness of the simulation 
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supercell, and this may cause some additional deviation 
as compared to DFPT. However, in our supercells the 
minimum q is quite smah (t^.03 bohr~^), and inspec- 
tion of the typical structured] of e°° as a function of q 
reveals that the deviations to be expected are in the order 
of 1 %. 

In closing, we discuss the case of non-polar materi- 
als: in such a case, the procedure outlined above does 
not apply, since no macroscopic polarization can occur 
in a system containing only a non-polar material. How- 
ever, since the internal fields in the superlattice layers 
are proportional to AP, we can set up a superlattice by 
alternating layers of the unpolarized material of interest 
and layers of some appropriate polarized material. The 
latter effectively function as a polarization supply for the 
unpolarized layer. In practice, to compute the dielectric 
constant of Si, we first calculate the dielectric constant 
and zero-field polarization for some polarized material, 
say SiC in the wurtzite structure. We then stack along 
the (111) direction of zincblende (i.e. (0001) of wurtzite), 
a superlattice such as [.../SiC/Si/SiC/Si/...]. Since SiC 
is polarized and we use periodic boundary conditions, the 
interface charge is 



(7) 



whence the dielectric constant of Si is trivially extracted. 
SiC should have the in-plane lattice constant of Si in or- 
der to avoid strain effects in the Si layer. Clearly its cal- 
culated properties in this specific, hypotetical realization 
are irrelevant: what counts is that it provides the polar- 
ization to create an interface charge and a depolarizing 
field inside Si. The above scheme, it turns out, is more 
of conceptual interest than of practical use for non-polar 
solids: interface states occur fairly easily at heterovalent 
junctions, spoiling the applicability of Eq. (0). Apply- 
ing the scheme in practice to non-polar solids will require 
quite some trial and error to identify a "clean" interface, 
and we did not pursue this further here. 

Summary and acknowledgements - In conclusion, our 
results indicate that the polarization-based method can 
produce theoretical dielectric constants within ^^5 to 10 
% of experiment, and is as accurate as DFPT. While 
much more limited in its general scope, our method ap- 
pears to be a useful alternative to DFPT for this kind of 
calculations. 

We thank Elena Manca for help in the perovskite cal- 
culations, and acknowledge special support from INFM 
within the PAISS program. 
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TABLE IL Electronic dielectric constants of several polar 
insulators. Our calculated values are compared with theoret- 
ical DFPT values and with experiment. 





Present 


DFPT 


Exp. 


GaAs 


12.53 


12.3" 


10.9-' 


AlAs 


9.37 


9.2" 


8.2-'^ 


InN 


8.49 




8.40^' 


w-SiC 


7.07 


7.28' 


6.65' 


GaN 


5.69 


5.41^^ 


5.70" 


ZnO 


4.65 




4.60-^ 


AIN 


4.61 


4.62^* 


4.68*" 


BeO 


3.15 




2.99' 


LiF 


2.19 




1.96^ 


PbTiOa 


8.28 


8.24"= 


8.64-' 


CaTiOs 


5.87 




5.8P 



" Ref. a ''Rcf. Ill '^Ref. Jl8|, '^Ref.Jl9|, '^Ref. |, ^Ref. |||, 
sRcf. ^ '^Ref. ^Ref. ||, ^Rcf. |J 



TABLE L Theoretical structural parameters for the mate- 
rials being studied. 





Structure 


ao (bohr) 


co/ao 


Uo 


GaAs 


zincblende 


10.60 






AlAs 


zincblende 


10.62 






LiF 


rocksalt 


7.50 






AIN 


wurtzite 


5.82 


1.619 


0.380 


GaN 


wurtzite 


6.04 


1.634 


0.376 


InN 


wurtzite 


6.66 


1.627 


0.377 


BeO 


wurtzite 


5.00 


1.610 


0.377 


ZnO 


wurtzite 


5.98 


1.616 


0.376 


PbTiOa 


cubic 


7.34 






CaTiOa 


cubic 


7.20 
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FIG. 1. Present theoretical £°° vs. experimental values. 
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FIG. 2. Relative percental error in e°° for DFPT and the 
present method (see text). 
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